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Abstract 

We obtain the reflection coefficients from a graphene sheet deposited on a material substrate 
nnder a condition that graphene is described by the hydrodynamic model. Using these coefficients, 
the gradient of the Casimir force in the configuration of recent experiment is calculated in the 
framework of the Lifshitz theory. It is shown that the hydrodynamic model is excluded by the 
measurement data at the 99% confidence level over a wide range of separations. From the fact 
that the same data are in a very good agreement with theoretical predictions of the Dirac model 
of graphene, the low-energy character of the Casimir interaction is confirmed. 
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I. INTRODUCTION 


Graphene is a two-dimensional sheet of carbon atoms which attracts much experimental 
and theoretical attention due to unique physical properties and great promise for various 
applications.- The most often used approach describing the electric and optical properties 
of graphene is the Dirac model.- It is applicable at low energies up to a few eV and assumes 
the linear dispersion relation for massless quasiparticles, which move with a Fermi velocity 
rather than with the speed of light.- The Dirac model was applied by many authors to calcu¬ 
late the van der Waals, Casimir and Casimir-Polder interactions in layer systems including 
graphene (see, for instance, Refs.-“-). The most straightforward formalism for performing 
these calculations uses the polarization tensor in (2-1-1 )-dimensional space-time.-^— Recently 
an equivalence of the formalisms exploiting the polarization tensor and the density-density 
correlation function has been proven.^® Furthermore, the formalism of the polarization ten¬ 
sor was appliedl^^i^ for comparison with the measurement data of recent experiment^ and 
demonstrated a very good agreement. 

Another approach used in theoretical description of the properties of graphene is the 
hydrodynamic model.— This model considers graphene as an inhnitesimally thin positively 
charged sheet, carrying a homogeneous fluid with some mass and negative charge densities. 
In the framework of this model, the dispersion relation for quasiparticles in graphene is 
quadratic with respect to the momentum. The hydrodynamic model was also considered 
and applied to calculate the van der Waals, Casimir and Casimir-Polder interactions in many 
papers.— It was foundi^ that in the interaction of graphene with either a Si or a An 
plate the hydrodynamic model predicts larger magnitudes of the Casimir free energy than 
the Dirac model. 

In this paper, we compare theoretical predictions of the hydrodynamic model with the 
experimental data of recent measurement^ of the gradient of the Casimir force between a 
Au-coated sphere and a graphene sheet deposited on a Si02 film covering a Si plate. For this 
purpose, we derive exact expressions for the reflection coefficients of the electromagnetic os¬ 
cillations on a three-layer structure, where one layer is a two-dimensional sheet described by 
the hydrodynamic model, whereas two other layers are described by the frequency-dependent 
dielectric permittivities. Then, the Casimir force and its gradient are calculated by using the 
standard Lifshitz theory^>^ in the proximity force approximation^ (PFA). We demonstrate 
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that theoretical predictions of the hydrodynamic model are excluded by the measurement 
data at the 99% conhdence level over a wide region of separations between the sphere and 
the graphene sheet. This allows to conclude that the hydrodynamic model of graphene does 
not describe such physical phenomena, as the van der Waals and Casimir forces. 

The paper is organized as follows. In Sec. II we derive the reflection coefficients for a 
graphene-coated substrate under an assumption that graphene is described by the hydro- 
dynamic model. Using these reflection coefficients, in Sec. Ill we calculate the gradient of 
the Casimir force in the experimental configuration of recent experiment^ and compare the 
theoretical results with the experimental data. Section IV contains our conclusions and 
discussion. 


II. REFLECTION COEFFICIENTS IN THE HYDRODYNAMIC MODEL 


We consider the amplitude reflection coefficients from a graphene sheet deposited 

on a thick plate (semispace) made of an ordinary material. Let us denote the reflection 
coefficient from a freestanding graphene sheet by and from a semispace in vacuum by 
Then, for the transverse magnetic (TM) and transverse electric (TE) polarizations of 
the electromagnetic field one obtainsl^ 
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Here, the signs minus and plus should be chosen for the TM and TE polarizations, respec¬ 
tively. 

In the framework of the hydrodynamic model, the reflection coefficients for a graphene 
sheet in vacuum calculated at the Matsubara frequencies along the imaginary frequency axis 
take the 
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Here, the Matsubara frequencies are = 27rkBTl/h, ks is the Boltzmann constant, T is 
the temperature, I = 0, 1, 2, qf = k‘]_ + k = |fcj_|, and fcj_ is the projection 

of the wave vector on the plane of graphene. The quantity K is the the single parameter 
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characterizing graphene in the framework of the hydrodynamic model. It has the meaning of 
the wave nnmber of a graphene sheet and is determined by the parameters of the hexagonal 
structnre of graphite (one 7r-electron per atom, resnlting in two vr-electrons per hexagonal 
cell). Calculation leads to^^— 
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( 3 ) 


where e and m are the charge and mass of vr-electrons and n = A/ (3\/3/^) with I = 1.421 A 
being the side length of the hexagon in a crystal lattice. The wave number in Eq. (j3|l 
corresponds to the frequency ojk = cK = 2.02 x lO^'^ rad/s. 

The reflection coefficients from the boundary plane of a semispace described by the 

dielectric permittivity Su = are the well known Fresnel coefficients 
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where 

ku = hi{i^i,k±) = kl + eii^. (5) 

Substituting Eqs. (|2]) and (jl]) in Eq. ([T]), we hnd the reffection coefficients from a graphene 
sheet deposited on a semispace made of ordinary material 
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For computational purposes, we express the reflection coefficients 
dimensionless variables 


( 6 ) 

in terms of the 
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where a is a parameter having the dimension of length (in the next section a has the meaning 
of separation distance between a graphene-coated substrate and a sphere). Then one arrives 
at 
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where 


kl = = j/2 + (£„ - l)Cf, K = 2aK. 


( 9 ) 


III. COMPARISON OF THE HYDRODYNAMIC MODEL WITH THE MEA¬ 
SUREMENT DATA 


In the first experiment on the Casimir effect in systems inclnding graphene, the gradient 
of the Casimir force was measured between a Au-coated hollow glass sphere of radius R = 
54.1 /im and a graphene sheet deposited on a Si02 him covering a Si plate.— The thickness 
of a Si02 him was D = 300 nm. The thickness of a Si plate (500 pm) was large enough to 
consider it as a Si semispace when calculating the Casimir force. In a similar way, the An 
coating on a sphere resulted in the same Casimir force as an all-An sphere. Measurements 
were performed by means of dynamic atomic microscope operated in the frequency-shift 
technique.—^— The force-distance relations were obtained with diherent applied voltages 
(20 repetitions) and with applied compensating voltages (22 repetitions) over the separation 
region from 224 to 500 nm for two diherent graphene samples. All the mean gradients of 
the Casimir force were found to be in a very good mutual agreement in the limits of the 
experimental errors.— As an example, in Fig. l(a,b) typical mean gradients of the Casimir 
force (the hrst sample, the measurement results obtained with applied compensating voltage) 
are shown as crosses at diherent separations a between the sphere and the plate. The vertical 
arms of the crosses indicate twice the total error AF' = 0.64 pN/m in measurements of the 
gradient of the Casimir force, and the horizontal arms are twice the error Aa = 0.4 nm in 
measurement of absolute separations. These errors were found at the 67% conhdence level, 
i.e., the true values of the force gradients and separations with a probability of 67% belong 
to the intervals [F'{a) — AF', F'{a) + AF'] and [a — Aa, a + Aa], respectively. 

Using the Lifshitz theory and the PFA, the gradient of the Casimir force between a An 
sphere and a graphene sheet deposited on a Si02 him covering a Si plate (semispace) is given 
by 
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Here, T = 300 K is the temperature at the laboratory, the prime on the summation sign 
multiplies by 1/2 the term with I = 0 , and the dimensionless variables y and Q are introduced 
in Eq. ([7]). The reflection coefficients from a An semispace (which replaces a sphere in the 
PFA) are given by Eq. (|1]). In terms of dimensionless variables, they take the form 


where, in accordance to Eq. 
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The reflection coefficients from a graphene sheet deposited on a Si 02 (fused silica) film 


covering a Si plate are expressed by the standard formulas of the Lifshitz theory between 
layered structures'*^ 
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Here, the reflection coefficients .Rtm te describe the reflection from a graphene sheet de¬ 
posited on a Si 02 semispace. If graphene is described by the hydrodynamic model, they are 
given by Eq. ([ 8 ]) with eu = /2a). The coefficients describe the reflection 

on the boundary plane between the two semispaces made of Si 02 and Si. They are the 
standard Fresnel reflection coefficients: 
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where 821 = e^^'\ic(li/2a) and k 2 i is defined similar to Eq. ([9]) with a replacement of eu with 
£ 21 - 

The quantity = e^^'^\ic(i/2a) entering Eq. flTTll is found^*^ using the Kramers- 
Kronig relation from the measured optical datai^i for extrapolated to zero frequency 

either by the Drude model with the plasma frequency ujp = 9.0 eV and relaxation parameter 
7 = 0.035 eV or by the nondissipative plasma model. Note that the above values of the Drude 
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parameters are in a very good agreement with the measured optical data.— Contrary to the 
expectations, the most precise experiments on measuring the Casimir interaction between 
metallic surfaces22 i“— are in agreement with the theoretical predictions using the plasma 
model extrapolation of the optical data and exclude the theoretical results using the Drude 
model extrapolation. Deep physical reasons why the plasma model extrapolation of the 
optical data is in agreement with the most precise measurements and the Drude model 
extrapolation is excluded by them remain unknown. Here, we perform all computations 
using both extrapolations for An. We find that for a sphere interacting with a graphene- 
coated substrate, where graphene is described by the hydrodynamic model, the difference 
arising from using different extrapolations is rather small. This allows to include it in the 
magnitude of the theoretical error like it was done in the case of metal-graphene interaction 
computed using the Dirac model of graphene.— 

In order to calculate the reflection coefficients o one also needs the values of 
and imaginary Matsubara frequencies. The B-doped Si plate used in the 

experiments^ had a resistivity between 0.001 and 0.005 D cm. This corresponds^S to a charge 
carrier density between 1.6 x 10^® and 7.8 x 10^®cm“^, i.e., well above the critical density 
at which the dielectric-to-metal phase transition accurs.^^ Then one obtains for the plasma 
frequency's the values between 5 x 10^^ and 11 x lO^'^ rad/s and for the relaxation parameter— 
7 1.1 X 10^^ rad/s. These Drude parameters were used to extrapolate the optical data^ 

for to zero frequency by means of either the Drude or the plasma model. Finally, the 

dielectric permittivity of Si at the imaginary Matsubara frequencies was found by means of 
the Kramers-Kronig relation like this was done previously in the literature.— Different types 
of extrapolation for a Si lead to only minor differences in the computed force gradients in 
the experimental conhguration. This is also taken into account in the theoretical error. For 
the dielectric permittivity of Si 02 an accurate analytic expression's has been used. 

The theoretical force gradients using the hydrodynamic model of graphene were computed 
by Eqs. (0), (HOD, (ED, (ED, and (ED- The computational results were corrected for the 
presence of surface roughness which contribution does not exceed 0.1% in this experiment.— 
The computed gradients of the Casimir force are shown as blue bands in Fig. l(a,b) over 
the entire measurement range. The uncertainty in the values of Up of Si and the differences 
between the predictions of the Drude and plasma model extrapolations of the optical data 
for An and Si determine the theoretical error, which is taken into account in the width of the 
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bands. As is seen in Fig. 1, the theoretical description of graphene using the hydrodynamic 
model is excluded by the data at the 67% conhdence level over the entire measurement 
range from 224 to 500 nm. With increasing conhdence level, the error bars of the mean force 
gradients and separation distances are also increasing. Thus, at the 95% conhdence level 
the maximum increase of the error bars is by a factor of two.— Note that if the errors are 
determined at the 95% or 99% conhdence levels the true values of the force gradients and 
separations belong to the wider intervals [F'(a) — AF', F'{a) + AF'] and [a —Aa, a +Aa] with 
probabilities of 95% and 99%, respectively. Taking this into account, it can be seen that 
over the range of separations from 224 to 450 nm theoretical predictions of the hydrodynamic 
model are excluded by the data at a higher, 95% conhdence level. If we further increase the 
conhdence level up to 99%, it is easily seen that theoretical predictions of the hydrodynamic 
model are still excluded, but this time over a more narrow separation range up to 360 nm. 

It is interesting also to compare the theoretical predictions of the hydrodynamic model 
with theoretical predictions of the Dirac model over a wider separation region from 220 nm 
to 1 pm. In Fig. 2 the gradients of the Casimir force in the conhguration of an experiment^ 
computed in this paper using the hydrodynamic model (the dashed line) and using the 
Dirac model^'^ (the solid line) are shown as functions of separation. As is seen in Fig. 2, 
the predictions of the hydrodynamic model remain to be larger than the experimentally 
consistent predictions of the Dirac model. The physical reason why the hydrodynamic model 
is not suitable for theoretical description of the Casimir force in layered systems including 
graphene may be in the linear dispersion relation inherent to graphene at low energies. This 
property makes a big difference between graphene and all types of ordinary dielectrics and 
metals. 

IV. CONCLUSIONS AND DISCUSSION 

In this paper, we have compared the measurement results for the gradient of the Casimir 
force between a Au-coated sphere and a graphene-coated substrate^^ with theoretical predic¬ 
tions of the hydrodynamic model of graphene used by many authors in previous literature. 
For this purpose, the reflection coefficients from the three-layer structure, where the hrst 
layer is graphene described by the hydrodynamic model and two other layers are described 
by the frequency-dependent dielectric permittivities, have been obtained. It was shown that 
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the hydrodynamic model of graphene is exclnded by the measnrement data over the entire 
measnrement range from 224 to 500 nm at the 67% conhdence level. Over a narrower sep¬ 
aration region from 224 to 360 nm an exclnsion of the hydrodynamic model by the data at 
even higher 99% conhdence level is demonstrated. 

The same experimental dataj^ was recently shown^ to be in a very good agreement 
with theoretical predictions nsing the Dirac model of graphene. Keeping in mind that the 
Dirac model is applicable at energies below a few eV, the results of this paper provide 
additional arguments in favor of the low-energy character of the Casimir interaction. In 
future it would be interesting to apply the hydrodynamic model for theoretical description 
of the rehectivity of graphene at higher energies, outside the applicability region of the Dirac 
model, and perform a comparison with respective experimental results. 
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FIG. 1: (Color online) The experimental data for the gradient of the Casimir force between a Au- 
coated sphere and a graphene sheet deposited on a Si02 film covering a Si plate (the first sample) 
are shown as crosses plotted at a 67% confidence level over different separation regions (a) from 224 
to 350 nm and (b) from 350 to 500 nm. The respective theoretical predictions of the hydrodynamic 
model of graphene are indicated as the blue bands. 
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FIG. 2: The gradients of the Casimir force between a Au-coated sphere and a graphene sheet 
deposited on a Si02 film covering a Si plate are computed using the hydrodynamic model (the 
dashed line) and the Dirac model (the solid line) as functions of separation. 
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